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Abstract 

We discuss various aspects of the multiparticle production processes at 
the LHC energy range with emphasis on the collective effects associated 
with appearance of the new scattering mode, which corresponds to the re- 
flective scattering and its impact on multiparticle production processes. 
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Introduction 



Nowadays the LHC is collecting data and providing experimental results at the 
largest world energy y/s = 7 TeV. Along with realization of its discovery poten- 
tial, the LHC experimental program renewed interest to the well known unsolved 
problems providing deepened insights into those issues. In this context the multi- 
particle production studies bring us a clue to the mechanisms of confinement and 
hadronization. Confinement of a color (i.e. the fact that an isolated color object 
has an infinite energy in the physical vacuum) is associated with collective, coher- 
ent interactions of quarks and gluons, and results in formation of the asymptotic 
states, which are the colorless, experimentally observable particles. The inelastic 
processes involve large number of particles in the final state. On the other side, 
the experimental measurements often reveal high degrees of a coherence in the 
relevant observables. No doubt these collective effects are very important for un- 
derstanding of the nonperturbative collision dynamics. Such collective effects are 
associated, in particular, with unitarity regulating the relative strength of elastic 
and inelastic processes and connecting the amplitudes of the various multipar- 
ticle production processes. Unitarity for the scattering matrix is formulated for 
the asymptotic colorless on-mass shell states. There is no universal, generally 
accepted method to implement full unitarity in high energy scattering. A related 
problem of the absorptive corrections and of their sign has a long history of discus- 
sion (cf. El and references therein). However, a choice of particular unitarization 
scheme is not just a matter of taste. 

Long time ago the arguments based on analytical properties of the scattering 
amplitude were put forward [2] in favor of the rational form of unitarization. It 
was shown that this form reproduced correct analytical properties of the scattering 
amplitude in the complex energy plane much easier compared to the exponential 
form. This specific form of unitarization [3] can be related to the confinement of 
color in QCD [ffl. 

Correct account for the unitarity is essential for the minimum bias multipar- 
ticle production processes, correlations, anisotropic flows and studies of phase 
transitions. 

The important point is that the region of the LHC energies is the one where 
the new, reflective scattering mode [5] can be observed. Such a mode naturally 
appears when energy grows and the rational form of unitarization being exploited 

. This mode can be revealed at the LHC directly measuring cr e i(s) and a to t{ s ) 
0. The reflective scattering is correlated with the self-damping of the inelastic 
channels and leads to the asymptotically dominating role of elastic scattering, i. 
e. a e i(s)/a tot (s) — y 1 at s — > oo while the both cross-sections, a e i(s) and a to t(s), 
tend to infinity in this limit. 

Of course, the obvious questions appear on compatibility of the reflective scat- 
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tering mode with dynamics of the multiparticle production, in particular, with the 
growth of the mean multiplicity in hadronic collisions with energy. Many mod- 
els and the experimental data suggest a power-ilke energy dependence of mean 
multiplicity^. 

In this review we discuss relation of the reflective scattering with confinement 
and hadronization, in particular, we consider the effects of the reflective scatter- 
ing in the multiparticle production processes at the LHC energies and apply a 
rational ([/-matrix) unitarization method j3]| to consider correlations in the multi- 
particle dynamics. We perform model calculations of global characteristics in the 
multiparticle production processes such as mean multiplicity, average transverse 
momentum, two-particle correlations, elliptic flow and demonstrate its quantita- 
tive and qualitative agreement with the first data obtained at the LHC. It appeared 
that reflective scattering leads to the prominent effects in the global observables 
and particle correlations. 

1 Rational form of unitarization and confinement of 
colored degrees of freedom 

Unitarity or conservation of probability, which can be written in terms of the scat- 
tering matrix as 

SS + = 1, (1) 

implies an existence of the two scattering modes at high energies - shadowing and 
reflection. Saturation of the unitarity condition for the scattering matrix in hadron 
collisions at small impact parameters b, takes place when scattering acquires re- 
flective nature, i.e. S(s,b)\b=o — > — 1 at s — > oo. Here S(s,b) is the elastic 
scattering S-matrix in the impact parameter representation. Approach to the full 
absorption in head-on collisions — the limit S(s, b) \ b=0 — >■ at s — > oo — does 
not follow from unitarity itself and is merely a result of the assumed saturation 
of the black disk limit. On the other hand, the reflective scattering is a natural 
interpretation of the unitarity saturation based on the optical concepts in high en- 
ergy hadron scattering. Such reflective scattering can be traced to the continuous 
increasing density of the scatterer with energy, i.e. when density goes beyond 
some critical value relevant for the black disk limit saturation, the scatterer starts 
to acquire a reflective ability. The concept itself is quite general, and results from 
the S'-matrix unitarity saturation related to the necessity to provide the total cross 
section growth at s — > oo. This picture predicts that the scattering amplitude at 
the LHC energies is beyond the black disk limit at small impact parameters. The 

'Discussions of power-iike energy dependence of mean hadronic multiplicity and list of refer- 
ences to the older papers can be found in ll9l [T0l[m 
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consequences for the total, elastic and inelastic cross-sections have been discussed 
in Id. 

Unitarity for the elastic scattering amplitude F(s, t) can be written in the form 

1mF(a, t) = H el (s, t) + H inel (s } t), (2) 

where H e i ine i(s,t) are the corresponding elastic and inelastic overlap functions 
introduced by Van Hove [[SJ. Physical meaning of each term in Eq. © is evident 
from the graphical representation in Fig. 1. The functions H e i :ine i(s, t) are related 




Figure 1: Eq. \2\in the graphical form. 

to the functions fr e j,mej( s > Yia tne Fourier-Bessel transforms, i.e. 

s f°° 
Jo 

The elastic and inelastic cross-sections can be obtained as follows: 



(3) 



0- e l,inel( S ) ~ ^H e i Me i(s,t = 0). (4) 

As it was already noted, the reflective scattering mode appears naturally in the 
U -matrix form of unitarization, where the 2 — > 2 scattering matrix element in the 
impact parameter representation is the following linear fractional transform: 

l+iU(s,b) 

S{s,b) = - — (5) 

1 — iu(s, b) 

U(s, b) is the generalized reaction matrix, which is considered to be an input dy- 
namical quantity. The relation © is one-to-one transform and easily invertible. 

Now, we would like to discuss relation of the rational unitarization form of 
scattering matrix with the confinement property of QCD. According to the con- 
finement, isolated colored objects cannot exist in the physical vacuum and there 
is no room for objects like quark-proton scattering amplitude, since isolated color 
object has an infinite energy in the physical vacuum. 

The important assumption in the derivation of the consequences of unitarity is 
the completeness of a set of the asymptotic states, those states include colorless 
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degrees of freedom or hadrons only. It can be considered as a questionable one in 
the QCD era. It might be reasonable to claim that a set of hadronic states does not 
provide a complete set of states and unitarity in the sense discussed above could 
be violated (cf. lfT2lO . It was stated in this paper that the Hilbert space, which 
corresponds to colorless hadron states and is constructed using vectors spanned 
on the physical vacuum, should, in principle, be extended. At the present time 
it is often underlined that the vacuum state is not unique; i.e. the colored current 
quarks and gluons are in fact the degrees of freedom related to a different vacuum. 
Thus, the vacuum state may not be considered anymore to be the state of the lowest 
energy (ground state). 

At the moment we would like to demonstrate that inclusion of the states cor- 
responding to the confined objects (e.g. colored current or constituent quarks) 
into a set of the asymptotic states would lead to a rational form of S -matrix unita- 
rization provided those states satisfy a certain constraint treated as a condition for 
confinement. We propose to manage these states similar to consideration of the 
states with indefinite metric in quantum electrodynamics as it was performed by 
N.N. Bogolyubov lfT3ll . To construct the S matrix (which operates in the physical 
subspace) let us consider state vectors |$_) at t — > — oo and |$ + ) at t — > +00 
each being the sum of the two vectors 

= \<p+) + K> 
|$_) = \(p_) + 

where \(p±) corresponds to the physical states and \co±) - to the confined states. 
So, we have that \<p±) = V±\&±) and \u±) = (1 — V±)\&±), where V± are the 
relevant projection operators relevant for the initial and final states, respectively. 

The scattering operator S (defined as |$ + ) = <S |$_)) should not, in principle, 
conserve probability and obey unitarity condition since it operates in the Hilbert 
space which includes subspace where confined objects with an undefined norm 
reside. Next, let us to impose condition on the asymptotic vectors \u±): 

+ \u+) = 0. 

It means that in- and out- vectors corresponding to the states of the confined ob- 
jects are just the mirror reflections of each other. Those reflections can be associ- 
ated with the impossibility for confined objects to propagate outside the hadron. 
In quantum mechanics such solution corresponds to standing or stationary wave 
where on average net propagation of energy is absent . Thus, the rational form of 
unitary scattering operator S 

s = v+s[i + (1 - v+s)}- 1 = (1 - u){\ + uy\ 
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in the physical subspace, i.e. \(p + ) = S\(f-)) can easily be obtained, since 

\<p + )=V + S(\v + ) + \u_)). 

We have started with scattering operator S in the Hilbert space which includes 
physical and unphysical subspaces and arrived to an unitary scattering operator 
S in the physical subspace. Crucial assumption was a constraint for the states of 
confined objects from the unphysical subspace + \lo + ) = 0, which we as- 
sume to be related to the confinement condition. The imposition of this constraint 
is an equivalent to the statement that the scattering matrix in the unphysical sub- 
space is identical with -1. Thus, one can say that if scattering matrix in the space 
spanned over confined states is identically equal to -1, then scattering matrix in 
the physical space is unitary. It is plausible therefore to assume that unitarity can 
be straightforwardly connected to confinement. 

Rational or [/-matrix form of unitarization was proposed long time ago lfT4l in 
the theory of radiation dumping. Self-damping of inelastic channels was consid- 
ered in [fT5ll and for the relativistic case such form of unitarization was obtained in 
0). It was demonstrated that it can be used for construction of a bridge between 
the physical states of hadrons and the states of the confined colored objects. 



2 Reflective scattering at the LHC energies 

In the rational form of unitarization the inelastic overlap function h ine i(s,b) is 
connected with the function U(s,b) by the relation 

hnU(s,b) 

h - i ^ b) = \i- iU ( S ,b)r (6) 

and the only condition to obey unitarity is ImU (s,b) > 0. The elastic overlap 
function is related to the function U(s, b) as follows 

The form of U(s, b) depends on particular model assumptions, but for qualitative 
purposes it is sufficient that it increases with energy in a power-like way and de- 
creases with the impact parameter like a linear exponent or Gaussian. To simplify 
the qualitative picture, we consider also the function U(s,b) as a pure imaginary. 
At sufficiently high energies (s > s R ), the two separate regions of impact param- 
eter distances can be anticipated, namely the outer region of peripheral collisions 
where the scattering has a typical absorptive origin, i.e. S(s, fe) | 6> j r(s) > and the 
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inner region of central collisions where the scattering has a combined reflective 
and absorptive origin, S(s, b)\ b<R ^ < 0. 

We discuss now the impact parameter profiles of elastic and inelastic overlap 
functions for the scattering picture with reflection. The well known absorptive 
picture of high energy scattering corresponds to the black disk limit at s — > oo. 
This limit has already been reached in the head-on proton-antiproton collisions 
at Tevatron, i.e. h e i(s,b = 0) ~ 0.25 lfT6ll . Thus, one can expect that in the 
framework of absorptive picture, black disk limit will be reached also at b ^ 
at higher energies and the profiles of h e i(s, b) and h ine i(s, b) are to be similar and 
have a form close to step function (Fig. 2). We consider now energy evolution 




b=R(s) b b=R(s) 



Figure 2: Typical picture of impact parameter profiles of the elastic and inelastic overlap 
function in the absorptive approach at asymptotic energies. 

of the elastic and inelastic overlap functions which includes reflective scattering. 
With conventional parametrization of the [/-matrix the inelastic overlap func- 
tion increases with energies at modest values of s. It reaches its maximum value 
hi ne i(s, 6 = 0) = 1/4 at some energy s = sr and beyond this energy the reflective 
scattering mode appears at small values of b. The region of energies and impact 
parameters corresponding to the reflective scattering mode is determined by the 
conditions h e i(s, b) > 1/4 and h ine i(s, b) < 1/4. The unitarity limit and black disk 
limit are the same for the inelastic overlap function, but these limits are different 
for the elastic overlap function. The quantitative analysis of the experimental data 
lfT71 gives the threshold value: ^/sr ~ 2 TeV. The function h ine i(s, b) becomes pe- 
ripheral when energy increases in the region s > s R . At such energies the inelastic 
overlap function reaches its maximum at b = R(s) where R(s) ~ In s. So, in the 
energy region, which lies beyond the transition threshold, there are two regions in 
impact parameter space: the central region of reflective scattering combined with 
absorptive scattering at b < R(s) and the peripheral region of pure absorptive 
scattering at b > R(s). Typical pattern of the impact parameter profiles for elastic 
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and inelastic overlap functions in the scattering picture with reflection at the LHC 
energies is depicted on Fig. 3. 




b=R(s) b b=R(s) b 



Figure 3: Typical qualitative picture of impact parameter profiles of the elastic and in- 
elastic overlap function in the reflective approach at the asymptotically high energies. 

It should be noted that at the values of energy s > s R the equation U(s, b) = 1 
has a solution in the physical region of impact parameter values, i.e. S(s, b) = 
at b = R(s). This line is shown in the s and b plane in Fig. 4 alongside with the 
regions where elastic S-matrix has positive and negative values. The dependence 



region I, no reflective scattering | g^g ^=Q 



i S(s,b) >o 




s=s r 



reflective scattering with probability |S(s,b)|' 



Figure 4: Regions of positive (absorptive scattering) and negative (absorptive and reflec- 
tive scattering) values of the function S(s, b) in the s and b plane. 



of S(s, b) on impact parameter b at fixed energies (in the region s > s R ) is depicted 
on Fig. 5. 
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S(s,b) 



+1 



Region of reflective scattering b<R(s) 




-1 



Figure 5: Qualitative impact parameter dependence of the function S(s, b)for the ener- 
gies s > SR. 



The probability of reflective scattering at b < R(s) and s > s R is determined 
by the magnitude of \S(s, b)\ 2 ; this probability is equal to zero at s < sr and b > 
R(s) (region I on Fig. 4). The behavior of R(s) is determined by the logarithmic 
dependence R(s) ~ j? In s. It is consistent with the analytical properties of the 
resulting elastics scattering amplitude in the complex t-plane and mass M can be 
related to the pion mass. 

Thus, at the energies s > s R the reflective scattering will mimic presence of 
repulsive core in hadron and meson interactions as well. The generic geometrical 
picture at fixed energy beyond the black disc limit is described as a scattering off 
the partially reflective and partially absorptive disk surrounded by the black ring 
which becomes gray at larger values of the impact parameter. The evolution with 
energy is characterized by increasing albedo due to the interrelated increase of 
reflection and decrease of absorption at small impact parameters. Asymptotically, 
picture of particle collisions with small impact parameters looks like collisions of 
hard spheres. 

3 Multiparticle production in the U -matrix approach 

To consider multiparticle production in the U -matrix approach it should be noted 
first that 




(8) 



n>3 
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where U n (s, b) is a Fourier-Bessel transform of the function 

Un(s,t) = I /]T^( 4 )(V gi -p a -p 6 )C/*( gi ,...., gn ;p' a ,^). (9) 
Til I Qjt) 

J j=l ™ i=l 

U n (qx,....,q n ;p a ,p b ). 

Here the functions U n (qi, q n ;p a ,Pb) and £4(<?i, q n ;Pa>,Pb>) correspond to 
the ununitarized (input or "Born") amplitudes of the process 

a + b — > 1 + .... + n, 

and the process with the same final state and the initial state with different mo- 
menta p' a and p' b 

a' + b' -)■ 1 + .... + n, 

respectively. They are the analogs of the elastic [/-matrix for the processes 2 — > n. 
It is important to note that the functions U n (s, t) are real ones (but not positively 
defined, contrary to the functions U n (s, b)). The impact parameter b is the variable 
which is conjugated to the variable t, where t — (p a — p' a ) 2 and the following 
sum rule is valid for the function I(s, b, q) 

J ^I(s, b, q) = (n)(s, b)JmU(s, b). (10) 

The sum in the right hand side of the Eq. © runs over all inelastic final states 
|n) which include diffractive as well as non-diffractive ones. Graphically, these 
relations are illustrated in Fig. 6. 

Im U 




n>3 



Figure 6: Unitarity for the function ImU. 

Then the inclusive cross-section of the process ab — >■ cX has the following 
form |[T8l [T9l , which is similar to the expression for the total inelastic cross- 
section: 
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where I(s, b, q) is the Fourier-Bessel transform of the functions which are defined 
similar to Eq. ® but with the fixed momentum q and energy E of the particle c 
in the final state (Fig. 7). 



I(s,t,q) 




n>3 



Figure 7: Unitarity for the function I(s, t, q). 

The impact parameter b is related to the impact parameters of the secondary 
particles by the relation EOl 

n 

b = ^Xibj, (12) 

i=l 

where Xi stands for the Feynman variable x of the i-ih. particle. 

The general relations listed above reflect the unitarity saturation only and 
do not include other interaction dynamics. To provide further insight into the 
mechanism of multiparticle production additional assumptions on the quark-gluon 
hadron structure and their interaction dynamics should be adopted. It is described 
in the following two sections 

4 Transient state of matter in hadron interactions 

We assume that the transient states of matter in hadron and nuclei collisions have 
the same nature and originate from nonperturbative sector of QCD, associated 
with the mechanism of spontaneous chiral symmetry breaking (xSB) in QCD 
[12111 . Due to this mechanism transition of current into constituent quarks takes 
place, the latter ones are the quasiparticles whose masses are comparable with 
a typical hadron mass scale. These constituent quarks interact via exchange of 
the Goldstone bosons which are collective excitations of the condensate and are 
represented by pions (cf. e.g. 11221 ). The general form of the effective Lagrangian 
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(£qcd — > ^e//) relevant for description of the non-perturbative phase of QCD 
includes the three terms [23] 



A// - £ X + &I + ^C7- 



Here £ x is responsible for the spontaneous chiral symmetry breaking and turns 
on first. To account for the constituent quark interaction and their confinement the 
terms Ci and Cc are introduced. Ci and Cc do not affect the internal structure of 
the constituent quarks. 

The picture of a hadrons consisting of the constituent quarks embedded into 
quark condensate implies that overlapping and interaction of the peripheral clouds 
occurs at the first stage of hadron interaction 



Figure 8: Schematic view of initial stage of the hadron interaction. 

Nonlinear field couplings transform then some part of their kinetic energy into 
internal energy related to a mass [|24ll231l . As a result we assume that the massive 
virtual quarks appear in the overlapping region and the effective field is generated 
in that way. This field is generated by QQ pairs and pions strongly interacting with 
quarks. In their turn pions themselves are bound states of the constituent quarks. 
At this stage the part of the effective Lagrangian Cc is turned off (it turns on again 
in the final stage of the reaction) and interaction is described by Cj. Its possible 
form was discussed in [|26l . The transient phase (effective field) generation time 



where At int is the total interaction time. This assumption on the almost instan- 
taneous generation of the effective field has obtained a support in the very short 
thermalization time revealed in the heavy-ion collisions at RHIC [|27ll28l . 

This picture assumes deconfinement at the initial stage of the hadron colli- 
sions and generation of the effective field common for both hadrons. Such ideas 
were used in the model [291 which has been applied to description of the elastic 




At 



At tp < An 
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scattering. Massive virtual quarks play a role of scatterers for the valence quarks 
in elastic scattering while their hadronization leads to production of the secondary 
particles. To estimate the number of such scatterers one could assume that certain 
part of hadron energy carried by the outer condensate clouds is being released 
in the overlap region to generate the massive quarks. Then this number can be 
estimated as 

iV( s , b) oc (1 - D h c } ® Df = N (s)D c (b), (13) 

m Q 

where mq - constituent quark mass, (kq) - average fraction of hadron energy car- 
ried by the constituent valence quarks. Here the function describes the conden- 
sate distribution inside the hadron h, and b is an impact parameter of the colliding 
hadrons. Thus, N(s, b) massive quarks appear in addition to N = + Uh 2 va- 
lence quarks. In elastic scattering those quarks are transient ones: they are trans- 
formed back into the condensates of the final hadrons. Calculation of the elastic 
scattering amplitude has been performed in 11291 . However, valence quarks can ex- 
cite a part of the cloud of the virtual massive quarks and those will subsequently 
fragment into the multiparticle final states. Such mechanism is responsible for 
the particle production and should lead to correlations between secondary parti- 
cles which originate from the particular quark cluster resulting from the valence 
quark excitation. In this impact parameter picture the strong forward-backward 
multiplicity correlations should be expected. 

Another mechanism contributing to the multiparticle production is a direct 
(i.e. not induced by interactions with valence quarks) hadronization of the massive 
quarks. 



5 Multiparticle production mechanism 

In sections 5 and 6 we consider correlations arising in the average multiplicity and 
transverse momentum behavior. 

Remarkably, the existence of the massive quark- antiquark matter at the stage 
preceding hadronization seems to be supported by the experimental data obtained 
at CERN SPS and RHIC (see ffl and references therein). 

Since the quarks are constituent, it is natural to expect direct proportionality 
between a secondary particles multiplicity and number of virtual massive quarks 
appeared in collision of the hadrons with given impact parameter: 

(n)(s,b) = a(n hl +n h2 )N (s)D F (b) + PN (s)D c (b), (14) 

with constant factors a and /3 and 

D F (b) = D Q ® D c , 
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where the function Dq(b) is the probability amplitude of the interaction of valence 
quark, which is in fact related to the quark matter distribution in this hadron- 
like object ll29~ll . The mean multiplicity (n)(s) can be calculated according to the 
formula 

, w x (n)(s,b)hi ne i(s,b)bdb 

= A f m • (15) 

J h inel (s,b)bdb 

It is evident that the peripheral profile of hi ne i(s, b) associated with reflective scat- 
tering suppresses the region of small impact parameters and the main contribution 
to the mean multiplicity is due to the region of b ~ R(s). 

To make explicit calculations we model the condensate distribution D c (b) 
and the impact parameter dependence of the probability amplitude Dq(b) by the 
exponential forms with the different radii. Then the mean multiplicity 

(n)(s,b) = aN (s)exp(-b/R F ) + PN (s)exp(-b/R c ). (16) 
The function U (s, b) is chosen as a product of the averaged quark amplitudes 

JV 

U(s,b) = H(f Q (s,b)). (17) 

0=1 

This factorization originates from an assumption of a quasi-independent nature of 
the valence quark scattering, iV is the total number of valence quarks in the col- 
liding hadrons. The 6-dependence of (/q) related to the quark formfactor Fq(q) 
has a simple form (fg) oc exp(— rnqb/E). Thus, the generalized reaction matrix 
(in a pure imaginary case) gets the form ||29ll 



U(s,b) = ig 



I — i JV 

m Q 



exp(-M6/0, (18) 



where M = ^2g =1 vtiq. At sufficiently high energies where increase of the total 
cross-sections is quite prominent we can neglect the energy independent term and 
rewrite the expression for U(s, b) as 

U(s,b) = ig {s/m Q ) N/2 exp(-Mb/0- (19) 

After calculation of the integrals (TT5T) we arrive to the power-like dependence 
of the mean multiplicity (n) (s) at high energies 

(n)(s) = as Sp +bs Sc , (20) 

where 
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Figure 9: Energy dependence of mean multiplicity, theoretical curve is given by the 
equation (n)(s) = as s (a = 2.328, S = 0.201). 

There are four free model parameters, a, and Rf, Rc, an d the freedom in 
their choice is translated to a, b and 5f, 5c- The value of £ = 2 is fixed from 
the data on angular distributions 11291 and for the mass of constituent quark the 
standard value mg = 0.35 GeV was taken. However, fit to experimental data 
on the mean multiplicity leads to approximate equality 5p — 5c and actually 
Eq. ((20b is reduced to the two-parametric power-like energy dependence of mean 
multiplicity 

(n) = as* 5 , 

which is in good agreement with the experimental data (Fig. 9). Equality 5p — 5c 
means that variation of the correlation strength with energy is weaker than the 
power dependence and could be, e.g. a logarithmic one. From the comparison 
with the data on mean multiplicity we obtain that 5 ~ 0.2, which corresponds to 
the effective masses, which are determined by the respective radii (M = 1/R), 
Mc — M F ~ 0.3mg, i.e. M F ~ M c — m^. The value of mean multiplicity 
expected at the LHC energy (y/s = 14 TeV) is about 1 10. 

Now we would like to make a remark on the mean multiplicity in the impact 
parameter representation. As it follows from the formulas of the section 2, the 
n-particle production cross-section cr n (s, b) 

U n (s,b) 



Then the probability 



P.(.,6) S -^ = ^!%. (22) 
(T ine i{s,b) lmU(s,b) 
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Thus, we observe cancellation of the unitarity corrections in the ratio of the 
cross-sections a n (s, b) and cr ine i(s, b). Therefore the mean multiplicity in the im- 
pact parameter representation 

(n}(s,b) = J2nP n (s,b) 

n 

is not affected by unitarity corrections. However, the above cancellation of unitar- 
ity corrections does not take place for the quantity (n)(s). The main contribution 
to (n)(s) originates from peripheral values of the impact parameter. Thus, typical 
inelastic event at the LHC energies is the event with a nonzero value of the initial 
impact parameter of collision. 

6 Rotation of transient matter and energy depen- 
dence of average transverse momentum 

We are going now to evaluate energy dependence of the average transverse mo- 
mentum of produced particles and propose a possible mechanism leading to this 
dependence. The geometrical picture of hadron collision discussed above implies 
that at high energies and non-zero impact parameters the constituent quarks pro- 
duced in overlap region carry large orbital angular momentum. It can be estimated 
as 

L(s,b)cxb^D c (b). (23) 

Due to supposed strong interaction between the quarks this orbital angular mo- 
mentum will lead to a coherent rotation of the quark system located in the overlap 
region as a whole. This rotation is similar to rotation of the liquid where strong 
correlations between particles momenta exist OTTl . In what follows we argue that 
this collective coherent rotation would lead to the energy dependence of the aver- 
age transverse momentum and can explain experimentally observed rising behav- 
ior of this quantity. It should be noted that discovery of the deconfined state of 
matter has been announced by four major experiments at RHIC [32 1. Despite the 
highest values of energy and density have been reached, a genuine quark-gluon 
plasma QGP (gas of the free current quarks and gluons) was not observed. The 
deconfined state reveals the properties of the perfect liquid, being a strongly in- 
teracting collective state and therefore it was labeled as sQGP. Using similarity 
between the hadronic and nuclear interactions we assumed that transient state in 
hadron interactions is also a liquid-like strongly interacting matter. As it was al- 
ready mentioned, the presence of large angular momentum in the overlap region 
will lead to coherent rotation of quark-pion liquid. Of course, there should be 
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experimental observations of this collective effect and one of them is the directed 
flow in hadron reactions, with fixed impact parameter discussed in OTTl . It is not 
impossible task to measure impact parameter of collision in hadron reactions with 
the help of the event multiplicity studies. But effects averaged over impact param- 
eter can be measured more easily using standard experimental techniques. So, it 
is natural to assume that the rotation of transient matter will affect average trans- 
verse momentum of the secondary hadrons in proton-proton collisions. Let for 
beginning do not take into account the other sources of the transverse momentum 
and temporally suppose that all average transverse momentum is a result of a co- 
herent rotation of transient liquid-like state. Then the following relation can be 
invoked 

(p T )(s,b) = KL(s,b), (24) 

where L(s, b) is given by Eq. (1231) and k is a constant which has dimension of 
inverse length. It is natural to relate it with inverse hadron radius, k ~ l/Rh- To 
calculate the average transverse momentum (pr(s)) the following relation with 
(Pt){s, b) will be used: 

J bdb{n}(s,b)h inel (s,b) 

where hi ne i(s, b) is the inelastic overlap function Now calculating the respective 
integrals we obtain the power-like dependence of the average transverse momen- 
tum (p T ) (s) at high energies 

(Pt)(s) = cs &c } (26) 

where 



>c 



1 (■> s 



1 - 



2 V m Q R Ci 

The value of £ is obtained from the data on angular distributions ||29l and for the 
mass of constituent quark the standard value vrtq = 0.35 GeV is taken. Of course, 
besides collective effects average transverse momentum would get contributions 
from other sources such as thermal distribution proposed long time ago by Hage- 
dorn [|33l . This part has no energy dependence and we take it into account by 
simple addition of the constant term to the power-dependent one, i.e.: 

(p T )(s) = a + cs Sc (27) 

Existing experimental data can be described well (cf. Fig. 10) using Eq. (1271) 
with parameters a = 0.337 GeV/c, c = 6.52 • 10" 3 GeV/c and 5 C = 0.207. 
The numerical value of Rc is determined by pion mass with better than 10% 
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Figure 10: Energy dependence of the average transverse momentum in pp-collisions, 
experimental data from 0?ll?5l/ . 

precision, R c — In the model the indices in the energy dependencies 

of average multiplicity and transverse momentum 5 and 5c are determined by 
the same expression and experimental data fitting with free parameters 5 and 5c 
confirms this coincidence with better than 10% precision also; note that the value 
5 = 0.201 follows from the experimental data analysis for the average multiplicity. 



7 Elliptic flow in proton collisions at the LHC ener- 
gies 



As it was already noted the most typical inelastic event at the LHC energies occurs 
with nonzero impact parameter. We consider therefore in this section peripheral 
hadronic collisions. The orientation of the reaction plane in pp-collisions can 
therefore be determined. 

There are several experimental probes of collective dynamics. Some of them 
were considered in the previous sections. Another observables are related to the 
anisotropic flows and among them the most widely discussed one is the elliptic 
flows 



which is the second Fourier moment of the azimuthal momentum distribution of 
the particles with a fixed value of px- The azimuthal angle <p is an angle of the 
detected particle with respect to the reaction plane, which is spanned by the colli- 



v 2 {pr) = (cos (20)) 



Pt ~ 




), 



(28) 
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sion axis z and the impact parameter vector b. The impact parameter vector b is 
directed along the x axis. Averaging is taken over a large number of the events. 
Elliptic flow can be expressed in covariant form in terms of the impact parameter 
and transverse momentum correlations as follows 

(b-p T ) 2 (b x p r ) 2 
v 2 (pt) = ( — -2 ) - ( -2 ), (29) 

FT FT 

where b = h/b. To get some hints on the possible behavior of the elliptic flow 
in proton collisions, it is useful to recollect what is known on this observable 
from nuclear collisions experiments. Integrated elliptic flow t> 2 at high energies is 
positive and increases with ^/snn- The differential elliptic flow ^(pt) increases 
with pt at small values of transverse momenta, then it becomes flatten in the 
region of the intermediate transverse momenta and decreases at large pr- It is also 
useful to apply a simple geometrical ideas which imply existence of the elliptic 
flow in hadronic reactions. Geometrical notions for description of multiparticle 
production in hadronic reactions were used by many authors, e.g by Chou and 
Yang in [36]. In the peripheral hadronic collisions the overlap region has different 
sizes along the x and y directions. According to the uncertainty principle we can 
estimate the value of p x as 1/ Ax and correspondingly p y ~ 1/ Ay where Ax and 
Ay characterize the size of the region where the particle originate from. Taking 
Ax ~ R x and Ay ~ R y , where R x and R y characterize the sizes of the almond- 
like overlap region in transverse plane, we can easily obtain proportionality of v 2 
in collisions with fixed initial impact parameter to the eccentricity of the overlap 
region, i.e. 

p2 _ d2 

The presence of correlations of impact parameter vector b and in hadron 
interactions follows also from the relation between impact parameters in the mul- 
tiparticle production: 

b = ^x,b,. (31) 

i 

Here Xi stand for Feynman xf of i-ih. particle, the impact parameters bj are con- 
jugated to the transverse momenta p^T- 

The above considerations are based on the uncertainty principle and angular 
momentum conservation, but they do not preclude the existence of the dynamical 
description, which will be discussed in the next section 

As it was noted in the section 6 an essential point of the proposed production 
mechanism is the rotation of the transient state due to the presence of a non-zero 
impact parameter in the collision. The unitarity saturation invokes a dynamical 
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mechanism leading to the peripheral nature of inelastic collisions at the LHC en- 
ergies. One can suggest therefore that the events with average and higher mul- 
tiplicity at the LHC energy y/s = 7 TeV correspond to the peripheral hadron 
collisions. Of course, the standard inclusive cross-section for unpolarized parti- 
cles being integrated over impact parameter b, does not depend on the azimuthal 
angle of the detected particle transverse momentum. 

When the impact parameter vector b and transverse momentum p T of the 
detected particle are fixed the function / does depend on the azimuthal angle 4> 
between vectors b and p T . The dependence on the azimuthal angle (f> can be 
written in explicit form through the Fourier series expansion 

1 x 
J(s,b,j/,pr) = —I (s,b,y,p T )[l + ^22v n (s,b,y,p T ) cosntfi]. (32) 

n=l 

The function I (s, b, £) satisfies to the following sum rule 

J I (s,b,y,p T )p T dp T dy = (n)(s,b)lmU(s,b), (33) 

where (n) (s, b) is the mean multiplicity depending on impact parameter. Thus, the 
bare anisotropic flow v n (s, b, y, p T ) is related to the measured flow v n as follows 

v n (s,b,y,p T ) = w(s,b)v n (s,b,y,p T ). 

where the function w(s,b) is 

w(s,b) = \l-iU{s,b)\- 2 . 

In the above formulas the variable y denotes rapidity, i.e. y = sinh _1 (p/m), 
where p is a longitudinal momentum. Thus, we can see that unitarity corrections 
are mostly important at small impact parameters, i.e. they modify anisotropic 
flows at small centralities, while peripheral collisions are almost not affected by 
unitarity. 

Now we are going to use the particle production mechanism described above 
for evaluation of the elliptic flow in pp-interactions. It was suggested that the ro- 
tation of transient matter will affect the average transverse momentum of the sec- 
ondary hadrons produced in proton-proton collisions. Going further, one should 
be more specific and note that, in fact, the rotation gives a contribution to the x- 
component of the transverse momentum and does not contribute to the y-component 
of the transverse momentum, i.e. 

Ap x = kL{s, b) (34) 
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while 

A Py = 0. (35) 

Assuming that p x = po + Ap x (Ap x <C p ) and p y = p , the contribution into 
the elliptic flow of the transient matter rotation can be calculated by analogy with 
the average transverse momentum calculation performed in |[3~71l . The resulting 
integral elliptic flow increases with energy 

Sr 

v 2 oc s , 

where Sc = 0.207. It should be noted here that transient matter consists of virtual 
constituent quarks strongly interacting via pion (Goldstone bosons) exchanges . 

The incoming constituent quark has a finite geometrical size determined by the 
radius tq and interaction radius Rq (Rq > Tq). The former one is determined by 
the chiral symmetry breaking mechanism and the latter one — by the confinement 
radius. Meanwhile, it is natural to suppose on the base of the uncertainty principle 
that size of the region where the virtual massive quark Q is knocked out from 
the cloud is determined by its transverse momentum, i.e. R ~ 1/pr- However, 
it is evident that R cannot be larger than the interaction radius of the valence 
constituent quark Rq. It is also clear that R should not be less than the geometrical 
size of the valence constituent quark tq for this mechanism be a working one. 
When R becomes less than tq, this constituent quark mechanism does not work 
anymore and one should expect vanishing collective effects in the relevant region 
of the transverse momentum. 

The value of the quark interaction radius was obtained under analysis of the 
elastic scattering [29] and it has the following dependence on its mass 

Rq = H m Q ~ l / m 7r (36) 

where £ ~ 2 and therefore Rq ~ 1 fm, while the geometrical radius of quark 
tq is about 0.2 fm. It should be noted the region, which is responsible for the 
small-pr hadron production, has large transverse dimension and the incoming 
constituent quark excites the rotating cloud of quarks with different values and di- 
rections of their momenta in that case. Effect of rotation will therefore be smeared 
off over the volume Vr and then one should expect that (Ap x )y R — 0. Thus, 

v%(p T ) = (v 2 ) Vr ~ (37) 

at small px- When we proceed to the region of higher values of px, the radius 
R is decreasing and the effect of rotation becomes more and more prominent, 
incoming valence quark excites now the region where most of the quarks move 
coherently, in the same direction, with approximately the same velocity. The mean 
value {Ap x ) Vn > and 

v^pr) ee (v^) Vn > (38) 
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and it increases with p T . The increase of v% with p T will disappear at R = Tq, i.e. 
at pt > 1/?"q, and saturation will take place. The value of transverse momentum 
where the flattening starts is about 1 GeV/c for tq ~ 0.2 fm. At very large 
transverse momenta the constituent quark picture would not be valid and elliptic 
flow vanishes as it was already mentioned. 

We discussed elliptic flow for the constituent quarks. Predictions for the el- 
liptic flow for the hadrons depends on the supposed mechanism of hadronization. 
For the intermediate values of pr the constituent quark coalescence mechanism 
[|38l l39l would be dominating one. In that case the hadron elliptic flow can be 
obtained from the constituent quark one by the replacement v 2 —> nyv® and 
Pt — > Pr/ n v^ where n v is the number of constituent quarks in the produced 
hadron. 

Typical qualitative dependence of elliptic flow in pp-collisions in this approach 
is presented in Fig.l 1 




P_T 



Figure 1 1 : Qualitative dependence of the elliptic flow v 2 on transverse momentum in 
pp- collisions. 

The centrality dependence of the elliptic flow is determined by the orbital 
angular momentum L dependence on the impact parameter, i.e. it should be de- 
creasing towards high and low centralities. Decrease toward high centralities is 
evident since no overlap of hadrons should occur at high enough impact parame- 
ters. Decrease of v 2 toward lower centralities is specific prediction of the proposed 
mechanism based on rotation since the central collisions with smaller impact pa- 
rameters would lead to slower rotation or its complete absence in the head-on 
collisions (it is also result of the symmetry in head-on collisions). Qualitative de- 
pendence of the elliptic flow on the impact parameter is similar to the one depicted 
in Fig.l 1 where variable b is to be used instead of transverse momentum. 
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8 Two-particle correlations in ^-collisions 



The ridge structure was observed first at RHIC in the two-particle correlation func- 
tion in the near-side jet production. It was demonstrated that the ridge particles 
have a narrow A(p correlation distribution (where <p is an azimuthal angle) and 
wide Ai] correlations (77 is a pseudorapidity). The ridge phenomenon was associ- 
ated with the collective effects of a medium li40ll . 

The similar structure in the two-particle correlation function was observed by 
the CMS Collaboration [41J. This is rather surprising result because the ridge 
structure was observed for the first time in pp-collisions. Those collisions are 
commonly treated as the "elementary" ones under the heavy-ion studies and there- 
fore often used as the reference process for detecting deconfined phase formation 
on the base of difference between pp- and AA-collisions. It is evident now that 
such approach should be revised in view of this new and unexpected experimental 
result. 

In the proposed explanation of the ridge effect at the LHC energy y/s = 7 TeV 
a dynamical selection of peripheral region in impact parameter space responsible 
for the inelastic processes is the important point. As it was already noted, the geo- 
metrical picture of hadron collision at non-zero impact parameters implies that the 
generated massive virtual quarks in overlap region could obtain very large initial 
orbital angular momentum at high energies. Due to strong interaction between the 
quarks this orbital angular momentum leads to the coherent rotation of the quark 
system located in the overlap region as a whole in the xz-plane. This rotation is 
similar to the rotation of liquid. The assumed particle production mechanism at 
moderate transverse momenta is an excitation of a part of the rotating transient 
state of massive constituent quarks (interacting by pion exchanges). Due to the 
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Figure 12: Interaction of the constituent quark with rotating quark-pion liquid. 
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fact that the transient matter is strongly interacting, the excited parts should be lo- 
cated closely to the periphery of the rotating transient state otherwise absorption 
would not allow to quarks and pions leave the interaction region (quenching). 

The mechanism is sensitive to the particular direction of rotation and to the 
rotation plane orientatation. This will lead to the narrow distribution of the two- 
particle correlations in A0. However, two-particle correlation could have broad 
distribution in polar angle (A77) in the above mechanism (Fig. 12). Quarks in 
the exited part of the cloud could have different values of the two components of 
the momentum (with its third component lying in the rotation a^-plane) since the 
exited region Vn has significant extension. 

Thus, the ridge-like structure observed in the high multiplicity events by the 
CMS Collaboration can be considered as experimental manifestation of the coher- 
ent rotation of the transient matter in hadron collisions. The narrow two-particle 
correlation distribution in the azimuthal angle is the distinctive feature of this 
mechanism. 

There should be other experimental observations of this collective effect. One 
was mentioned already. It is the directed flow v\ in hadron reactions, with fixed 
impact parameter discussed in OTTl . Rotation of transient matter will affect also 
elliptic flow i> 2 and average transverse momentum of the secondary particles pro- 
duced in proton-proton collisions. 

The ridge effect has triggered flow of possible explanation using different 
mechanisms ES|43IM|45l|4a|47ra 

The one described in this paper and in Il53ll relates appearance of the ridge 
effect with the reflective scattering at the LHC energies. Absence of the reflective 
scattering at lower energies explains why the ridge effect was not observed at the 
ISR or Tevatron. This important point on the ridge energy dependence needs to be 
explained by any particular mechanism proposed as a source of the two-particle 
correlations. 

9 Genuine QGP formation 

The dynamical mechanism of the average transverse momentum growth, anisotropic 
flows v\, t> 2 , two-particle correlations (ridge) originate from the collective effect 
of transient matter rotation, while dynamics of average multiplicity growth is re- 
lated to the mechanism where a nonlinear field couplings transform the kinetic 
energy to internal energy. Formation of a genuine quark-gluon plasma in transient 
state in the form of the noninteracting gas of free quarks and gluons would result 
in disappearance of the effects of rotation of the transient state. Vanishing the en- 
ergy dependent contribution to the average transverse momentum should therefore 
be expected, i.e. average transverse momentum would reach maximum at some 
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energy and then will decrease till eventual flat behavior (cf. Fig. 13). Vanishing 




Collison energy 



Figure 13: Qualitative energy dependence of the average transverse momentum in pp- 
collisions in case of genuine QGP formation at super high energies. 

rotation will lead to vanishing anisotropic flows v\ and t> 2 . Of course, this pic- 
ture is a rather qualitative one since the nature of phase transition from strongly 
interacting matter (associated with quark-pion liquid in the model) to the genuine 
quark-gluon plasma is not known. Currently, there is no experimental indications 
that such a transition will indeed take place at super high energies. At the same 
time there are theoretical arguments that existence of genuine QGP (ideal gas) 
would contradict to the confinement prop[erty of QCD. This contradiction arises 
in the many dynamical mechanisms of hadronization Il54l[55tt . However, the issue 
of the genuine QGP existence can only be resolved by the experimental searches. 

It should be noted that inverse phase transition (from parton gas in parton 
model to liquid) could explain a saturation phenomena in deep inelastic processes 

m. 

The natural question appears about the role of orbital angular momentum in 
those mentioned above phase transitions. 

10 Spin correlations due to reflective scattering in 
the case of genuine QGP formation 

It is interesting to make a conclusion on the possible effect of the imbalance of the 
orbital angular momentum in the case of the genuine QGP formation. It should 
be noted that we can consider separately particles production in the forward and 
backward hemispheres II571 . Let us consider for example particles produced in 
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the forward hemisphere. The orbital angular momentum in the initial state can be 
estimated as 

U = (39) 
2 2 

where i?(s) is the interaction radius. The orbital angular momentum in the final 
state is then 

/ 4 2' 

where we have taken into account that (n)(s)/2 particles with the average frac- 
tion of their longitudinal momentum (xl) (s) are produced at the impact parameter 
R(s)/2 due to reflective scattering. The average fraction of longitudinal momen- 
tum (xl)(s) according to the hypothesis of limiting fragmentation Il58ll would not 
decrease with energy. Thus we arrive to the negative imbalance of the orbital 
angular momentum 

AL = U - L f * - <">(«)(«*>(«)), (41) 

/ 2 2 v 2 h v ' 

2 2 2' 
This negative AL should be compensated by the total positive spin S of final 
particles (since the particles in the initial state are unpolarized) 

S = -AL (43) 

This spin alignment of produced particles appears in the reflective scattering mode 
when the particles are produced in the region of impact parameters b = R(s). 

The vector of spin S lies in the transverse plane but it cannot be detected 
through the transverse polarization of single particle due to azimuthal symmetry 
of the production process (integration over azimuthal angle (p). However, this 
effect can be traced measuring transverse spin correlations of two particles (siSj). 
The most evident way to reveal this effect is to perform the measurements of 
the spin correlations of hyperons whose polarizations can be extracted from the 
angular distributions of their weak decay products. 

Spin correlation should be stronger for the light particles and the weakening 
is expected for heavy particles since they should be produced at smaller values of 
impact parameters. 

Thus, we can expect appearance at the LHC energies of strong spin correla- 
tions of final particles as a result of the prominent reflective scattering mode in 
the case if genuine QGP is formed in the transient state. Vanishing anisotropic 
flows, decreasing average transverse momentum and appearance (simultaneous) 
of the secondary particles polarization are thus the signals of the genuine QGP 
formation. 
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11 Effects of the reflective scattering mode for nu- 
clear collisions 



In this final section we proceed from proton collisions to the collisions of nuclei. 
Consider central collision of two identical nuclei having N nucleons in total with 
center of mass energy y/s per nucleon and calculate nucleon density ur(T, /j) = 
N/V in the initial state at given temperature T and baryochemical potential [i 
in the presence of the reflective scattering. The effect of the reflective scattering 
of hadrons is equivalent to decrease of the volume of the available space which 
the hadrons are able to occupy in the case when reflective scattering is absent. 
Thus followings to van der Waals method, we must then replace volume V by 
V ~ Pr{s)Vr{s)^, i.e. we should write 

V - p R {s)V R (s)f 

where n(T, //) is hadron density without account for reflective scattering and 
Pr(s) is the averaged over volume Vr(s) probability of reflective scattering: 

^00 = 7tVt [ \S(s,r)\ 2 d 3 x. 



V R (s) 



Vr(s) 



The volume Vr(s) is determined by the radius of the reflective scattering. Here 
we assume spherical symmetry of hadron interactions, i.e. we replace impact 
parameter b by r and approximate the volume V R (s) by V R (s) ~ (47r/3)i? 3 (s). 
Hence, the density n R (T,ix) is connected with corresponding density in the ap- 
proach without reflective scattering n(T, /j) by the following relation 

irp \ n(T,/i) 
nR{T,ii) - 



1 + a(s)n(T, jj) 



where a(s) = pr(s) Vr(s)/2. Let us now estimate change of the function ur(T, /i) 
due to the presence of reflective scattering. We can approximate Pr(s) by the 
value of \S(s, b = 0)| 2 which tends to unity at s — > oo. It should be noted that 
the value ^/sr ~ 2 TeV IfTTll . Below this energy there is no reflective scatter- 
ing, a(s) = at s < sr, and therefore corrections to the hadron density are 
absent. Those corrections are small when the energy is not too much higher than 
sr. At s > Sr the value of a(s) is positive, and presence of reflective scattering 
diminishes hadron density. We should expect that this effect would already be 
noticeable at the LHC energy y/s ~ 5 TeV in Pb + Pb collisions. At very high 
energies (s — > oo) 

n R (T,fj,) ~ l/a(s) ~ M 3 /ln 3 s. 
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Figure 14: Phases of strongly interacting matter. 



This limiting dependence for the hadron density appears due to the presence of 
the reflective scattering which results in similarity of head-on hadron collisions 
with scattering of hard spheres. It should be noted that this dependence has been 
obtained under assumption on spherical symmetry of hadron interaction region. 
Otherwise, limiting dependence of the hadron density in transverse plane can only 
be obtained, i.e. transverse plane density of hadrons would have then the following 
behavior 

tir{T, y) ~ M 2 / In 2 s. 

Thus, the following chain of transitions between different (physical and non- 
perturbative) vacuum states can be foreseen as the temperature increases at the 
constant value of chemical potential fi: 

|0) p / l (Hadron gas) — > |0) np (Quark-pion liquid) — > |0) p / l (Hadron liquid). 

The corresponding phase diagram depicted in Fig. 14. 

Thus, the reflective scattering mode can help in searches of the deconfined 
state and studies of transition mechanism to this state of matter. 
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